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Abstract 

We study the CP-violating terms of the spin-spin correlations in neutralino pro- 
duction and their subsequent two-body decays into sleptons plus leptons at the 
ILC. We analyze CP-sensitive observables with the help of T-odd products of 
the spin-spin terms. These terms depend on the polarizations of both neutrali- 
nos, with one polarization perpendicular to the production plane. We present 
a detailed numerical study of the CP-sensitive observables, cross sections, and 
neutralino branching ratios in the Minimal Supersymmetric Standard Model 
with complex parameters. 
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1 Introduction 



It has been pointed out that the amount of CP violation in the Standard Model (SM) is 
not sufficient to explain the baryon-antibaryon asymmetry of the universe [1], and that 
additional sources of CP violation are required [2] . Many extensions of the SM can give 
rise to such sources of CP violation. The violation of the CP symmetry is an interesting 
topic in its own right and deserves a diligent consideration. Supersymmetric (SUSY) 
extensions of the SM provide new sources of CP violation, as they include several 
new parameters which can be complex. For instance, in the neutralino sector of the 
Minimal Supersymmetric Standard Model (MSSM) two complex parameters appear, 
which lead to CP-violating effects in reactions involving neutralinos. These parameters 
are the higgsino mass parameter \i = \/i\e 1 ^, and the U(l) gaugino mass parameter 
Mi = \M 1 \e l ' /)1 , given in the usual parametrization of modulus and phase. 

These phases, on the other hand, contribute to the electric dipole moments (EDMs) 
of electron, neutron, and that of the atoms 199 Hg and 205 T1 [3], and it is found in 
general that for phases of the size 0(1), the EDMs are beyond their experimental 
upper bounds. However, the extent to which the EDMs can constrain the CP phases 
also depends on most of the other model parameters, and thus strongly depends on 
the considered model, see e.g. Refs. [3,4]. 

In this respect the high-luminosity e + e~ International Linear Collider (ILC) is 
considered an ideal machine to perform precision measurements, in order to determine 
the model parameters of the MSSM with the required accuracy [5]. In neutralino 
production and decay at the ILC, it has been shown which CP-even observables are well 
suited to access the CP- violating MSSM parameters [6,7]. However to directly prove 
CP violation in the MSSM, and to determine the CP-violating phases unambiguously, 
a measurement of CP-odd observables is obligatory. 

In this paper, we study CP-sensitive observables in neutralino production, 

e + + e-^x-+xh i,j = 2,3,4, i^j, (1) 

based on T-odd correlations [8] which appear in the spin-spin correlation terms of the 
amplitude squared. These terms involve the polarizations of both neutralinos, with 
one polarization perpendicular to the production plane. Such a normal polarization 
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component is a genuine signal of CP violation (neglecting higher order effects due to 
final state interactions [8]). The polarizations of the neutralinos can be analyzed in 
their decays, that's why we consider the leptonic channels 

X? + y° 3 -i't R + t' T , £,£> = e,». (2) 

Due to angular momentum conservation, the decay distributions of the final leptons E, 
£' are correlated to each other, and thus allow us to probe the spin-spin correlations. 

In a previous publication, we have analyzed in this way the CP-sensitive spin-spin 
correlations for chargino production and decay [10]. Other works done on CP-sensitive 
observables in neutralino pair production at the ILC have taken into account the 
decay of only one neutralino, where again the normal polarization component signals 
CP violation [11, 12]. Even the potential of transverse beam polarizations for CP 
observables in neutralino production has been analyzed [13,14]. CP observables have 
also been studied in decays of neutralinos, which originate from sfermions [15]. 

The paper is organized as follows. In Section [2} we define the Lagrangians and 
complex couplings for neutralino production. In Section [31 we present the analytical 
formulae for the amplitude squared of neutralino production and decay. In Section HI 
we identify the T-odd products in the spin-spin terms of the amplitude squared. In 
Section [5], we define the CP-sensitive observables which probe these terms. We present 
numerical results in Section EJ where we also estimate the measurability of the CP- 
sensitive observables. We give a summary and the conclusions in Section [7J 



2 Lagrangians and complex couplings 

In the MSSM, neutralino production e + e~ —>■ XiX^ proceeds via Z boson exchange in 
the s-channel, and selectron cl,r exchange in the t- and w-channels, see the Feynman 
diagrams in Fig. [TJ The Lagrangians for production and decay are [6, 16] 

Czee = g —Z^[L e P L + R e P R ]e } (3) 

cos v w 

1 Note that generally parity-conserving neutralino decays, like x° — ► Z\\ or %° — > fix®, would lead 
to vanishing CP-sensitive observables. Due to the Majorana properties of the neutralinos, the left 
and right neutralino couplings to the Z (and Higgs) have equal absolute values, and thus all spin- 
and spin-spin correlations would be lost, see the discussion in Ref. [9]. 
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3 Cross section 



The differential cross section for neutralino production e + e — > an d decay x° — * 
~£%r^i X°j -» ^'l,r^, can be written 

da = — |T| 2 dLips , (11) 

with the center-of-mass energy y/s, and the Lorentz invariant phase space element 
dLips, see Appendix The amplitude squared |T| 2 was calculated in Ref. [6] in the 
spin density matrix formalism^! 



4|A(x°)| 2 |A(^ 



P A Dj + J2 £p ^% D j 



a=l 



6=1 a,b=l 



(12) 



with the neutralino propagators A(y. ,) = l/[» 2 — rn 2 + im v o T v o 1. The ampli- 

v l)J/ ' 1 XiJ XiJ X iJ X i,j' 

tude squared has contributions from neutralino production (P) and decay (D). The 
terms P and D i: Dj are those parts of the spin density production and decay matrices, 
respectively, that are independent of the polarizations of the neutralinos. The con- 
tributions Sp and H%. depend on the polarization basis vectors s" of neutralino xl- 
Similarly, S P and H b D depend on the polarization basis vectors of neutralino 
See Appendix |Bj Eq. ( IB. 10[) for the explicit definition of the spin vectors. We choose 
a coordinate frame such that a, b = 3 denote the longitudinal polarizations, a, b = 1 
the transversal polarizations in the production plane, and a, b = 2 the polarizations 
normal to the production plane. The decay terms D±, Dj, S^., and H b D are given in 
Appendix [A] The full expressions for the production terms P, S p , S P , and £ P can 
be found in Ref. [6]. 

The contributions to the amplitude squared which depend on the polarizations of 
both neutralinos are the spin-spin correlation terms £ P . The CP-sensitive parts of 
the spin-spin correlation terms include one neutralino spin vector with a component 
perpendicular to the production plane, i.e., ab = 12, 21, 23, 32 [6] 

Xf(ZZ) = ^^\A(Z)\\Rlc R + Llc L )lm{0^0':f*} f ab , (13) 



2 For a detailed discussion of the spin density matrix formalism, we refer to Ref. [18]. 
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Xf(Ze L ) = J L e c L A(Z)[A* u (e L ) + A* t (e L )} M/dtfW> f ab , (14) 
2 cos^ Vw 

Z P b (e L e L ) = -^c L A u (e L )A* t (e L ) Im{(^) 2 (/^) 2 } f ab , (15) 

T${Ze R ) = - 9 " R e c R A(Z)[At(e R ) + A* t (e R )) \^{f^0^} f ab , (16) 
2 cos z tiw j j 

Z a P b (e R e R ) = £c R A u (e R )A* t (e R ) Im{(/*) 2 (/5*) 2 } f ab . (17) 

The dependence on the longitudinal beam polarizations is given by 

c L = (l-V-)(l + V+), c R = (l + V-)(l-V + ) , (18) 

with V- and V+ the degrees of longitudinal polarization of the electron and positron 
beam, respectively, with — 1 < V± < 1. Generally the contributions from the exchange 
of e R (e^) are enhanced and those of (e R ) are suppressed for V- > 0, V+ < (V- < 
0,V + > 0). The propagators are A(Z) = i/(s — m 2 z ), A t (e LjR ) = i/(t — m\ L ), 
A u (e L . R ) = i/(u-ml L J, with s = (p e - + p e +) 2 , t = (p e - -p Xj ) 2 , u = (p e - -p Xi f [6]. 

The spin-spin correlation terms £p in Eqs. (U3l) - ffTTj) explicitly depend on the imag- 
inary parts of the products of neutralino couplings, Im{0^ L 0-^*}, Im{/^/^*0^ 1 '}, 
Im{(^) 2 (/ e ^) 2 }, lm{f e R J e fO':f}, and Im{(/«) 2 (/^) 2 }. For i ? j they are mani- 
festly CP-sensitive, i.e., sensitive to the phases 0^ and <pi of the neutralino sector. 
These imaginary parts of the couplings are multiplied by T-odd factors f ab , which we 
discuss in detail in the next section. 

4 T-odd products of the spin-spin correlations 

The kinematical dependence of the spin-spin correlation terms of neutralino produc- 
tion, Eqs. f|T3l) - f|T7j) . is given by the T-odd function [6] 

f ab = (Pe+-P Xj )[Pe-,P X ^S^S h J + (p e ^p Xi )[p e+ ,p X]} S^S h X] ] 

+ (Pe+-S Xj )[Pe-,P Xt ,P Xj ,S* i } + (p e --slJ{p e+ ,p Xt ,p X] ,S b J , (19) 

with the short hand notation of the epsilon product of the four four-vectors pi 

\PuP2,P3,Pa] = S^afi Pi P\ P3 Pi > With £ 0123 = ~1- (20) 
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Since each of the spacial components of the four-momenta or spins changes sign under 
a naive time transformation, t — > — t, the epsilon product, and thus the function f ab , 
is T-odd. In Appendix [Bj we give f ab also in the laboratory system. 

In order to identify the T-odd products which appear in the spin-spin correlations 
of production and decay, we analyze the corresponding terms of the amplitude squared, 
Eq. (fl2i) . in more detail 

\T\ 2 D £p ^ oc £ f ab ■ (p r s a Xi ) ■ (pe-s b x .) , (21) 

a,b=l a, 6=1 

where the scalar products (pe ■ s£.) and (p? ■ s^.) stem from Efj. and E^,., respectively, 
see Eqs. (IA.2j) or flA.4j) in Appendix [S] Using the explicit expression for f ab , Eq. ffl9l) . 
and the completeness relation for the neutralino spin vectors, Eq. (IB. lip , the right- 
hand side of of the second equation in Eq. f[2~Tj) can be written as 

Ot = (Pe+-P Xj )\Pe-,P X vP£iPe'] + (Pe--P X i)\Pe+,P X jiPeiP£'] 

+ (Pe+-Pt')\Pe-,PxnPxjiPt] +(Pe--P£)\Pe+,PxiiPxiiPt'] ■ ( 22 ) 

We have now identified the CP-sensitive terms of the neutralino spin-spin correlations. 
They are proportional to the T-odd product Ot, Eq. (1221) . which can now be used 
to define various CP asymmetries and CP observables in neutralino production and 
decay. Due to their similar kinematical dependence, the definition of CP observables 
is analogous to those in chargino production and decay, see Ref. [10]. 



5 CP-sensitive observables 

In this Section, we define various CP-sensitive observables, which depend on the T- 
odd parts of the spin-spin correlations for neutralino production and decay. For an 
operator O, we define its expectation value by [10] 

/0|T| 2 dLips If da 

If the operator O is chosen of the form like the T-odd terms Ot, Eq. (122]) . the CP- 
sensitive parts of the spin-spin correlations in neutralino production can be projected 
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out 

/ O Ef E a D Z b D dLips 
<°> = S PdX dLips ' < 24 > 

with an implicit sum over (a, b) = (1,2), (2, 1), (2, 3), (3, 2). In the numerator remain 
only the CP-sensitive parts of the spin-spin terms of the amplitude squared. Only they 
contain the T-odd product O. In the denominator, all spin- and spin-spin correlation 
terms vanish, and only the spin-independent part P Di Dj contributes. Note that for 
the phase space element dLips in Eq. (I!MI) . we have already used the narrow width 
approximation for the propagators, see Eq. (]C.22j) . 



In general the largest observables are obtained by using an operator O, which 
exactly matches the kinematical dependence of the CP-sensitive terms in the amplitude 
squared, that is O = Ot, Eq. (1221) . In the literature, this technique is sometimes 
referred to optimal observables [19]. Thus for the operator Ot we define the two 
CP-sensitive observables 

(O t ) and A T = (Sgn(0 T )> . (25) 

Neglecting higher order effects due to final state interactions [8], the observable At is 
a CP asymmetry. It is the expectation value for the sign of the T-odd product Ot, 
and can be written as 

At = N +~ N - , (26) 
N+ + N- ' K J 

the difference of the number of events with positive (N + ) and negative (iV_) sign of 
Ot, normalized by the total number of events N = N + + AL. On the other hand, 
(Ot) is the expectation value of the momentum configuration Ot itself for the event 
sample. 

Two further T-odd products were considered in Ref. [10]. One product is ob- 
tained from Ot, Eq. ( |22i) . in replacing the four-momenta by the (normalized) three- 
momentum vectors in the center-of-mass system, see Appendix (Bj 

Ot = (Pe- • Pt') Pe~ ■ (P Xl X Pi) + (Pe- ' Pi) Pe~ ' (P X , X Pi') , ( 27 ) 



S 



with p = p/\p\. In contrast to Ot, Eq. (1221) . this product does not involve the energies 
of the neutralinos and leptons. For the operator Ot, we again define two CP-sensitive 
observables 



(O t ) and A T = (Sgn(O r )) . (28) 

Since both T-odd products Ot and Ot include the neutralino momentum p Xi and/or 
p Xj , their experimental reconstruction is required. For the subsequent two-body decays 
of the neutralinos which we consider here, the neutralino momentum three-vectors can 
be reconstructed up to a sign ambiguity in their second component, if the masses of 
the involved particles are known, see for example Ref. [13]. 

A T-odd product which does not depend on the neutralino momenta is obtained 
by substituting on the right hand side of Eq. ( 1271) the neutralino three-momenta by 
the corresponding decay lepton three-momenta p Xi — > pe and p x , — > p# [10], 

0't = Pe- ■ (Pe + Pf) p e - ■ (Pe x Pe) ■ (29) 
Also for 0' T we define a CP-sensitive observable and its corresponding asymmetry, 

(0' T ) and A' T = (Sgn(d' T )) . (30) 



Thus, depending on the type of correlation used, two classes of CP observables can 
be defined; one class requires the reconstruction of the neutralino momenta, the other 
class not. However, as we will show in the numerical section, the largest observables 
are obtained if indeed the neutralino momenta can be reconstructed. 



5.1 Relative signs of the CP observables 

Each of the above defined CP observables depends in principle on the various decay 
channels of the two neutralinos. For each neutralino, these are 

xl - ft + r, (si) 

- h+P, ( 32 ) 
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for £ = e, /i, and also 



t£ + t 



(33) 




(34) 



if the decay into the (usually) heavier left slepton is kinematically allowed. However, 
only the sign of the CP observable changes, depending on the charge and the type (L or 
R) of the two decay sleptons, for an overview see Table HJ The reason is that the signs 
of the corresponding two neutralino decay terms, E 1 ^ and T1 D ., only depend on the 
charge and the type of the two decay sleptons, see Eqs. flA.2j) and (1A.4I) . respectively 
The absolute value of an observable is independent of the particular decay channels, 
since the absolute values of the couplings \fg' i ' R \ or \ff- R \ of the decay sleptons, as well 
as their masses, cancel in the numerator and denominator of Eq. (1241) . This means 
in turn that we have to distinguish from which neutralino Xi or X°j the final state 
leptons £ and £' originate. Without that information, the contributions to the CP 
observables from the final leptons with charge combinations £~£'~ and £ + £' + would 
identically cancel those contribution from £ + £'~ and £~£' + . 

Furthermore if also the decay into £l is kinematically possible, the type of the 
sleptons, £l or £ R , into which the neutralinos decay, has to be determined. Such a 
discrimination can in principle be accomplished by using the different energy distribu- 
tions of the decay leptons, since their kinematical limits depend on the mass difference 
of the decaying neutralino and slepton. 

Note however, that if the final lepton momenta are assigned correctly, one is able 
to reconstruct the production plane. Although the two lightest neutralinos in the 
end of the decay chains, \°k Xi^ carry away their missing momentum, 

the ambiguities in the azimuthal angles of the produced neutralinos can be resolved 
with a measurement and correct assignment of the four final lepton momenta, see the 
discussion in Ref. [13]. Certainly the feasibility of such an event reconstruction can 
only be answered by a detailed experimental analysis, which is however beyond the 
scope of the present work. 
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Table 1: Relative signs of the CP-sensitive observables for different decay combinations 
of neutralino x° — * ^l(r)^ T (t°P row )> an d neutralino x° _ » ^l(_r)^' T (^ e ^ column). 
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+ 




it 


+ 
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5.2 Theoretical statistical significances 

We have defined various kinds of CP-sensitive observables, which are based on the dif- 
ferent T-odd products O = Tl Tl 0' T . They either include (0 T) O t ) or not include 
(O f T ) the neutralino momentum. In order to be able to compare these observables 
quantitatively, we define their theoretical statistical significances. A comparison of 
the numerical values of (O) and A = (Sgn(0)) alone cannot be used to decide which 
observable is more sensitive to the CP phases. In addition, we are sometimes facing 
situations where large CP observables and asymmetries correspond to processes with 
small neutralino production cross sections or branching ratios, and vice versa. Such 
effects can be considered by combining both the CP observable and the cross section 
into one statistical quantity. 

We define the theoretical statistical significance of the CP observable (6), where 
6 = O, or 6 = Sgn(C), by [13,20] 

S\0] = VN • (35) 

For neutralino production and decay the number of events is 

N = F N x C x a(e+e- - x°X?) x [BR(x° - e~+ e~) x BR(*° - e+ e~) 

+BR(x° - e+ e") x BR(x? -> e\ e~) 
+BR(x° - e+ e") x BR(x° - e£ e") 
+BR(x° - e~+ e") x BR(x° - e~+ e")] , (36) 
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with the integrated luminosity C. The combinatorial factor F N takes into account 
all possible neutralino decays into sleptons with different flavors and charges. We 
assume that the branching ratios of the neutralinos do not depend on them, i.e., 
BR(x£ - e+ e") = BR(x° k - e~ e+) = BR(x° k - K \T) = BR(x° k - K /i+), for 
n = L and R. The combinatorial factor is thus F N — 4 x 4 = 16, if we sum the two 
lepton flavors e, /i and the two charges, an d -C- 

The statistical significance S is equal to the number of standard deviations to which 
the corresponding CP observable can be determined to be non-zero over statistical fluc- 
tuations. For example, 5 = 1 implies a measurement at the statistical 68% confidence 
level, assuming an ideal detector, and full reconstruction of signal and background. 
Thus our definition of S is theoretically motivated, and can only be regarded as as an 
upper bound on the confidence level which at best can be obtained. In order to give 
realistic values of the statistical significances, a detailed experimental study would be 
required, which is however beyond the scope of the present work. 

Also higher order corrections have to be included in a comprehensive analysis. 
Although we expect the influence of electroweak corrections to our observables and 
asymmetries to be small, the corrections to the neutralino masses and production cross 
sections can be 10% at one-loop level [22]. The neutralino branching ratios for two- 
body decays may receive CP-even one-loop corrections of up to 16% in some cases [23]. 
For chargino production additional CP-sensitive terms contribute at higher order to 
the production cross section, which has recently been discussed in Ref. [24]. 



6 Numerical results 

We present numerical results for the CP-sensitive observables and asymmetries for 
neutralino production e + e~ — > X2X3, an d decay xt ~^ ' ■> X3 ~^ ^r^' t -, f° r — 
e, /J. We will study the dependence of the CP observables on the phases and moduli 
of the higgsino and U(l) gaugino mass parameters \x = |/z|e^ M and M± = jMije^ 1 , 
respectively, in the framework of the general MSSM. In this model the restrictions 
on the phases from the electron and neutron EDMs are less severe compared to the 
constrained MSSM [4]. Thus we will not take the EDMs into account, and show the 
full phase dependence of the observables. 
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Table 2: Benchmark scenario for e + e — > X2X3) an d decay %° ~ * ^r^ T ) X3 ~^ ^'r^i 
for £, £' = e,fi, at i/i = 500 GeV with beam polarizations (V-,V+) = (0.9, —0.6). 



M 2 




M 0i tan /3 


m 


270 GeV 


200 GeV 


0.5vr 3 


80 GeV 



The results are calculated with a center-of-mass energy of y/s = 500 GeV. We 
choose longitudinal beam polarizations (V-,V + ) = (0.9, —0.6), which enhance the e# 
exchange contribution. The feasibility of measuring the observables also depends on 
the neutralino production cross section and decay branching ratios, which we discuss in 
detail. For a comparison of the CP observables, and for giving an upper bound on the 
confidence levels, we also present a theoretical estimate of their statistical significances. 

For the neutralino widths and branching ratios, we include the two-body decays [12] 

X° - L + i, Oi + vi, Z + x° m , h + x° m , W ± + xt, (37) 

with m < i; k — 1, 2; n = R, L for I = e, fi, and n = 1, 2 for £ = t. We neglect three- 
body decays. We use the GUT inspired relation \M%\ = 5/3M 2 tan 2 9w, such that the 
dependence of the CP observables on the modulus of Mi is investigated by using M 2 . 
In order to reduce the number of free parameters further, we parametrize the slepton 
masses by M 2 , and m = 80 GeV fixed, which enter in the approximate solutions 
to the renormalization group equations, see Appendix |A] We take stau mixing into 
account, and fix the trilinear scalar coupling parameter A T = 250 GeV. We use the SM 
parameters sin 2 6 W = 0.2315, m w = 80.41 GeV, m z = 91.187 GeV, a = 7.8125 x 10~ 3 . 

The CP-sensitive neutralino coupling factors in Eqs. (113j) - (117p are zero for i = j, or 
vanishing phases M and 0i. They are largest for 0i = 0.57T (or l.bir), and for a strong 
gaugino-higgsino mixing M 2 w \fi\. We find that a small value of tan/? is preferred to 
have large CP observables and large significances. Therefore we center our numerical 
discussion around a scenario with tan j3 = 3 and a strong gaugino-higgsino mixing. The 
parameters are summarized in Tab. [21 The corresponding particle masses, branching 
ratios, and the cross section are listed in Tab. [3j For this scenario, we analyze the 
phase dependence of the CP observables, and then their dependence on and M 2 . 
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Table 3: SUSY masses, neutralino branching ratios and production cross section, for 
the benchmark scenario. The branching ratios are summed over £ — e,fi and both 
slepton charges. 



m v o = 121 GeV 

Al 


m e - 8 = 157 GeV 


BR(x° -> 4£) = 66% 


m v o = 171 GeV 


m SL = 256 GeV 


BR(x° -> f x r) = 34% 


m v o = 207 GeV 


m fl = 157 GeV 


BR(x° -> = 66% 


m ± = 158 GeV 

Al 


m f2 = 256 GeV 


BR(x° 3 - nr) = 34% 


m + = 318 GeV 

-A2 


ms = 246 GeV 


v{e + e- - xlxl) = 79 fb 



6.1 Phase dependence 



In Fig. 2(a) , we show the <j>x dependence of the CP asymmetries At ( |25l) . At ( 1281) . and 
A' T (1301) . The asymmetries vanish at <pi = 0, ir, 2ir, where the neutralino couplings are 
real. They obtain largest values at 0i ~ 0.57T and 4>i ~ 1-57T of about ^4 r = ±19%, 



At = ±16%, and A' T = ±8%. In Fig. 2(b), we show the same asymmetries as a func- 



tion of 0^, setting <pi = ir. They show a similar behavior, and again At attains the 
largest values of all three asymmetries. We do not present plots of the corresponding 
observables, since they show similar phase dependences as their corresponding asym- 
metries. They obtain values of (O t ) = -2.56 x 10 11 GeV 6 , (6 T ) = -0.062, and 
(p' T ) = 0.027, for the scenario defined in Tab. EJ 

In order to compare now the CP asymmetries A with their corresponding CP ob- 
servables (O), we show their theoretical significances S as a function of 0i in Fig. [3j 



First we observe that the observables, Fig. 3(b) , have generally larger significances than 



their counterpart asymmetries, Fig. 3(a). Secondly, the observables and asymmetries 
which are based on the T-odd products Ot and Ot, which include the neutralino 
momentum, have the largest significances. They would be best suited for measuring 
CP phases in the neutralino spin-spin correlations. The significance of (Ot) is twice 
as large as that of (0' T ). However, for their measurement a reconstruction of the 
neutralino momenta, i.e., the production plane is necessary, which will be experimen- 
tally more involved. The need to only reconstruct the final state leptons might be an 
advantage in a realistic experimental environment. However, a detailed investigation 
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Figure 2: Dependence of the CP asymmetries At (dotted), At (dashed), and 
A' T (solid), (a) on the phase 0i with 0^ = 0, and (b) on the phase M with <f>\ = ir, 
and the other parameters as defined in Tab. EJ 

S[0] 

30: 

25 ] 
20: 




Figure 3: Phase- dependence of the significances of (a) the asymmetries At (dotted), 
At (dashed), A' T (solid), and (b) of the observables (Ot) (dotted), (Ot) (dashed), 
(O't) (solid), with an integrated luminosity of £ = 500 fb -1 , for the scenario as defined 
in Tab. [2j 

which observable will be best suited can only be answered by a thorough experimental 
analysis, which is beyond the scope of the present work. In order to further illus- 
trate the different magnitudes of the asymmetries At and A T , we show them and the 
corresponding significances as a function of the phases M and 0i in Fig. HI 

Finally, it should be noted that a measurement of observables which depend only 
on cj)^ will be helpful to disentangle CP-violating effects in the neutralino system, 
which could originate both from 4>i and M . This could be possible by investigating 
CP observables in the chargino system [10,21] which solely depend on M . Finally a 
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Figure 4: Contour lines of the asymmetries At and .A^, and their statistical signif- 
icances in the plane, with an integrated luminosity of C = 500 fb _1 , for the 
scenario as defined in Tab. [2J 
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global fit of CP-even [6,7] and CP-odd [11-14] observables in the neutralino system 
could allow for a complete determination of the phases. 



6.2 fi and Mi dependence 

In order to estimate the significances of the CP-sensitive observables in a larger region 
of the parameter space, we now analyze the neutralino cross sections, branching ratios 
and, as an example, the asymmetry At (|25l) in the |/i|-M 2 plane. 



In Figs. 5(a) and 5(b) , we show the neutralino branching ratios which are summed 
over both lepton flavors £ — e,// and charges, i.e., BR(x° — > £r£) = 4 x BR(x° — > 
e^e - ), for k — 2,3. In the gray shaded area, the chargino mass is m x ± < 100 GeV, 
and thus near or below the exclusion limit of LEP2 [25] . In region A, the neutralinos 
are below the decay threshold, rn x o^ < rn-^ and thus the corresponding two-body 
decays are closed. The neutralino x\ * s always lighter than £l in the shown region of 
the |//|— M% plane. We find that the x\ branching ratio into left sleptons is smaller than 



BR(%3 — > £ L £) < 1%. In Fig. 5(a) and 5(b) , the decay channels into the lightest Higgs 



and Z bosons open to the right of the dashed lines, which indicate the kinematical limit 
m x o a = m x o + mz, respectively. However, these channels would lead to vanishing CP 
observables, due to the Majorana properties of the Higgs and Z boson couplings to the 



neutralinos, as discussed in the introduction. Along the dotted contour in Fig. 5(b) , the 
decay channel xl ~^ W^xT opens, which also considerably reduces BR(xg — > £r£) to 
the right of that contour, for \fi\ > M 2 . The neutralino xX an d X3 branching ratios into 
staus become larger than those into selectrons for |//| > M 2 . If the tau momenta can be 
reconstructed, these decay channels can also be used to measure the CP observables. 
However due to stau mixing, the observables will be reduced compared to the decays 
into selectrons or smuons, see the discussion in Ref. [9]. 

The neutralino production cross section (T23 = a(e + e~ — > X2X3) is shown in 



Fig. 5(c), It reaches values up to 130 fb for M 2 ~ 250 GeV and w 150 GeV. 
In the region B, the neutralinos are too heavy and above the production threshold, 
m x o + m x o > y/s = 500 GeV. The combined cross section of production and decay, 



a = (T 2 3 x BR(x° — >■ £r£) x BR(x° — > £r£), is shown in Fig. 5(d). One can see the 
combination of the kinematically excluded regions from production and decay. The 
cross section a reaches up to 65 fb. 
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Figure 5: and M 2 dependence of (a) the neutralino branching ratio BR(%2 — > £r£), 
(b) the branching ratio BR(x|| — > £r£), (c) the neutralino production cross section 
<7 2 3 = cr(e + e~ — > X2X3)) an d (d) the combined cross section of production and decay, 
o" = 023 x BR(%2 —> £r£) x BR(%3 — > £r£), for the scenario as defined in Tab. [2j In 
region A the neutralinos are below the decay threshold, m x Q 23 < m i R i an d in region B 
they are above the production threshold, m x o + m x o > y/s = 500 GeV. In the gray 
shaded areas the chargino mass is m ± < 100 GeV. The dashed contours in (a), (b) 
indicate the kinematical limit ^ x « 3 = m x \ + mz > respectively- The dotted contour in 
(b) indicates the limit m x o = mw + m x * ■ 
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Figure 6: Contour lines of (a) the asymmetry At and (b) its statistical significance 
^.At] in the \fi\-M2 plane, for the scenario as defined in Tab. [3 In region A neutralino 
X2 is below the decay threshold, m x o < , and in region B the neutralinos are above 
the production threshold, m x o + m x o > y/s = 500 GeV. In the gray shaded areas the 
chargino mass is m + < 100 GeV. 



In Fig. [HI we show the asymmetry At and its corresponding significance 5'[^4r] in 
the \n\-M2 plane. The asymmetry reaches values up to —30%, while the significance 
goes up to 50 standard deviations close to the kinematical limit m x o = mj> R , at M2 ~ 
300 GeV and m 180 GeV. At that point, the asymmetry of the correlation 0' T , 
that does not need the reconstruction of the neutralino momenta, reaches A' T = 13%, 
which corresponds to a significance of about if>[y4y] = 25. 
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7 Summary and conclusions 



We have analyzed CP observables in neutralino production, which are sensitive to 
the physical phases of the gaugino parameter M±, and the higgsino parameter \i. 
The observables and asymmetries rely on T-odd products in the neutralino spin-spin 
correlations, which appear on tree-level. The CP-sensitive spin-spin correlations are 
those terms of the matrix element, which include the polarizations of both neutralinos, 
with one component normal to the production plane. These spin-spin correlations of 
the neutralinos can be analyzed via angular distributions of the decay leptons x° — > £ £. 

In order to probe the CP-sensitive spin-spin correlation terms, we have defined 
different T-odd products. One class only involves the final lepton momenta, which 
has the advantage that it is not necessary to reconstruct the production plane. The 
second class of T-odd products also includes the neutralino momenta. Based on these 
T-odd products, we have studied two sorts of CP-sensitive observables. One sort are 
CP-sensitive observables, which are the expectation values of the T-odd products. The 
other sort are their corresponding asymmetries, which give the expectation value of 
the sign of the T-odd products. 

In our numerical analysis for X2X3 production, we have found that the observables 
are largest in mixed scenarios with small tan/5. We have defined theoretical signif- 
icances to decide, which CP observable is most sensitive to the CP phases. For a 
linear collider with ^fs = 500 GeV and longitudinally polarized beams, (V-,V+) = 
(0.9,-0.6), with an integrated luminosity of C = 500 fb -1 , the CP-sensitive observ- 
ables that only include the momenta of the decay leptons yield theoretical significances 
of S < 25 for 0! = 0.57T. We find larger theoretical significances up to S < 50 for the 
CP-sensitive observables that need a reconstruction of the neutralino momenta. How- 
ever, only a detailed experimental study with background and detector simulations can 
show whether the CP-sensitive observables are accessible. We hope that our results 
motivate such a study. 
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Appendix 



A Decay terms D and E c 



D 



The coefficients in Eq. ([12]) of the neutralino decay matrices for the decay into right 
sleptons xt ~* ^r with £ — e, //, are [12] 

D k = yl/^Ko-mfJ, (A.l) 
^ = ^9Vi\ 2 m x o(s c Xk -p e ), (A.2) 

where the sign in parenthesis holds for the charge conjugated process xt ~~ > ^r ^~ • 
For the decay into the left sleptons xt ~ > ^ = e, A<, the coefficients are 

D k = 9 ^\fk\ 2 Ho-ml), (A.3) 
Z c Dk = M9 2 \fk\ 2 rn x o(s Xk ■ p t ) , (A.4) 

where the sign in parenthesis holds for the charge conjugated process Xk ~^ ^~ • 

In order to reduce the free MSSM parameters, we parametrize the slepton masses 
with an approximate solution to the renormalization group equations (RGE) [26] 

mj = m 2 + m 2 e + 0.23M 2 - m 2 z cos 2 f3 sin 2 9 w , (A.5) 
m | £ = m 2 Q + m 2 e +0.79M 2 + m 2 z cos2(3(~ + sm 2 6 w ) , (A.6) 
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m 2 ^ = m 2 + m] + 0.79M 2 2 + cos 2(3 , (A.7) 



with mo the common scalar mass parameter at the GUT scale. 



B Momentum and polarization vectors 

We choose a coordinate system with the z-axis along the p e - direction in the center- 
of-mass system. The four-momenta of the neutralinos x° and x° are 

p Xi . = g (^, T sin0,O, T cos0) , (B.8) 



with their energies and common momentum 

E ^ ~ — vi — ' q ~ — Vi — ' (R9) 



respectively, and the kinematic function A(a, b, c) = a 2 + b 2 + c 2 — 2(a6 + ac + 6c). The 
scattering angle is 9<t(p e - ,p Xj ), whereas the azimuthal angle can be set to zero, due 
to rotational invariance around the beam axis [27]. 

The three spin basis vectors of x°i and X°j are chosen to be 

\ ' Xi,j Xi,j ' ' 

4, = fo,^^) =(0,0,1,0), 

S L = Z^-U^r-^oj) =§^f^-,Tsin^,0, T cos^ . (B.10) 



They fulfill the orthonormality relations s£ • = — <5 C(i , s£ fe • p Xk = 0, and the 
completeness relation [6, 18] 



E°X-°Z = -r + ^- (b.ii) 



The four-momenta of the leptons in the decays x° — ► Z A and x° - ► ^' are 



P£ = (1, cos ^ sin ^, sin ^ sin ^, cos ^) , (B.12) 
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Pe — cos0£' sin^/, sin0^ sin^/, cos^/) , (B.13) 

respectively, with 

m 2 — m~ ml— m\ 

2{E Xi + q cos -de) 2(£ Xj - q cos V t < ) 

and the decay angles 

cos $e = sin 6 sin Oe cos 4>e + cos 9 cos 9i , 

cos'di' = sin 6 sin Off cos <figt + cos 6 1 cos Off . (B.15) 

With these definitions, the T-odd products f ab ( TT9i) of the spin-spin correlation terms 
in the laboratory system are 

f U = - l -E Xi sqsm 2 6 , f 21 = ±E Xj sq sin 2 , (B.16) 
/ 23 = ^m x .sgsin(20) , f 2 = ~m Xi sq sin(20) . (B.17) 



C Phase space 



The Lorentz invariant phase space element in Eq. ffTTj) is given by [12,27] 

dLipS = (2^p dLi P s ( s ' Pxi > Px J ) ds Xi dLi P s ( s Xi > PI) pi)ds Xi dLips(s x . , p^,, p^), 



(C.18) 



with s Xi ^ = p Xj . The different factors of the phase space element are 



1 q 

dLips(s, p Xi ,p Xj ) = -sm6 dO , (C.19) 

o7T v/ 5 
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dLi PS (s Xi , P/ , Pl) = 2(27r) 2 m 2 _ m 2 Sil1 °l d ^ d ^ > ( C - 20 ) 



|2 



dLip S (s Xi , p-^Pff) = 2(27r) 2 m 2^ m 2 Sil1 ^' d ^' d ^' • ( C - 21 ) 

23 



We use the narrow width approximation for the propagators in Eq. (|12p . 



which is justified for T/m <C 1, which holds in our case with T < O(lGeV). Note, 
however, that the naive (9(T/m)-expectation of the error can easily receive large off- 
shell corrections of an order of magnitude and more, in particular at threshold, or due 
to interferences with other resonant or non-resonant processes. For recent discussions 
of these issues, see Ref. [28]. 
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